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Abstract 

The production of } /xp mesons is studied in pp collisions at = 7 TeV with 
the CMS experiment at the LHC. The measurement is based on a dimuon sam- 
ple corresponding to an integrated luminosity of 314 nb -1 . The J / ip differential 
cross section is determined, as a function of the } /xp transverse momentum, in 
three rapidity ranges. A fit to the decay length distribution is used to separate 
the prompt from the non-prompt (b hadron to } / xp) component. Integrated over 
J/ip transverse momentum from 6.5 to 30 GeV/c and over rapidity in the range 
\y\ < 2.4, the measured cross sections, times the dimuon decay branching fraction, 
are 70.9 ± 2.1(stat.) ± 3.0(syst.) ± 7.8 (luminosity) nb for prompt // xp mesons assum- 
ing unpolarized production and 26.0 ± 1.4(stat.) ± 1.6(syst.) ± 2.9 (luminosity) nb for 
/ 1 xp mesons from b-hadron decays. 
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1 Introduction 

Heavy-flavour and quarkonium production at hadron colliders provides an important test of 
the theory of Quantum Chromodynamics (QCD). The production of ]/tp mesons occurs in three 
ways: prompt ]/ip produced directly in the proton-proton collision, prompt ]/ip produced in- 
directly (via decay of heavier charmonium states such as Xc), and non-prompt ]/ip from the 
decay of a b hadron. This paper presents the first measurement of the differential inclusive, 
prompt and non-prompt (b hadron) ]/ip production cross sections in pp collisions at a centre- 
of-mass energy of 7 TeV, in the rapidity range \y\ < 2.4, by the Compact Muon Solenoid (CMS) 
experiment. 

Despite considerable progress in recent years IQ3-0, quarkonium production remains puzzling 
and none of the existing theoretical models satisfactorily describes the prompt J/t/> differen- 
tial cross section [3-5 1 and polarization (6) measured at the Tevatron [7] . Measurements at the 
Large Hadron Collider (LHC) will contribute to the clarification of the quarkonium produc- 
tion mechanisms by providing differential cross sections in wider rapidity ranges and up to 
higher transverse momenta than was previously possible, and with corresponding measure- 
ments of quarkonium polarization. Cross-section results are largely dependent on the ]/xp po- 
larization, as different polarizations cause different muon momentum spectra in the laboratory 
frame. Given the sizeable extent of this effect, for prompt ]/ip mesons (where the polarization 
is presently not well described by the theoretical models) we choose to quote final results for 
different polarization scenarios, instead of treating this effect as a source of systematic uncer- 
tainty. 

Non-prompt ]/xp production can be directly related to b-hadron production, leading to a mea- 
surement of the b-hadron cross section in pp collisions. Past discrepancies between the Teva- 
tron results (both from inclusive |5| and exclusive [8[ measurements) and the next-to-leading- 
order (NLO) QCD theoretical calculations, were recently resolved using the fixed-order next- 
to-leading-log (FONLL) approach and updated measurements of the b — > ]/ip fragmentation 
and decay 0[lOl. Measured cross-section values and spectra are also found to be in agreement 
with Monte Carlo generators following this approach, such as MC@NLO MTTllT2f . 

The paper is organized as follows. Section 2 describes the CMS detector. Section 3 presents 
the data collection, the event trigger and selection, the }/ip reconstruction, and the Monte Carlo 
simulation. Section 4 is devoted to the evaluation of the detector acceptance and efficiencies 
to detect ]/rp events in CMS. In Section 5 the measurement of the ]/tp inclusive cross section is 
reported. In Section 6 the fraction of ]/ip events from b-hadron decays is derived, and cross- 
section results are presented both for prompt ]/xp production and for ]/ip production from b- 
hadron decays. Section 7 presents comparisons between the measurements and model calcula- 
tions. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diam- 
eter, providing a field of 3.8 T. Within the field volume are the silicon pixel and strip tracker, 
the crystal electromagnetic calorimeter and the brass/ scintillator hadron calorimeter. Muons 
are detected by three types of gas-ionization detectors embedded in the steel return yoke: Drift 
Tubes (DT), Cathode Strip Chambers (CSC), and Resistive Plate Chambers (RPC). The mea- 
surement covers the pseudorapidity window \rj\ < 2.4, where tj = — ln[tan(0/2)] and the 
polar angle 9 is measured from the z-axis, which points along the counterclockwise beam di- 
rection. The silicon tracker is composed of pixel detectors (three barrel layers and two forward 
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3 Data sample and event reconstruction 



disks on each side of the detector, made of 66 million 100 x 150 /mi 2 pixels) followed by mi- 
crostrip detectors (ten barrel layers plus three inner disks and nine forward disks on each side 
of the detector, with 10 million strips of pitch between 80 and 184 }im). Thanks to the strong 
magnetic field and the high granularity of the silicon tracker, the transverse momentum, pj, of 
the muons matched to reconstructed tracks is measured with a resolution of about 1 % for the 
typical muons used in this analysis. The silicon tracker also provides the primary vertex posi- 
tion, with ~ 20 }im accuracy. The first level (LI) of the CMS trigger system, composed of custom 
hardware processors, uses information from the calorimeters and muon detectors to select the 
most interesting events. The High Level Trigger (HLT) further decreases the rate before data 
storage. A much more detailed description of the CMS detector can be found elsewhere [13 j. 

3 Data sample and event reconstruction 

3.1 Event selection 

The analysis is based on a data sample recorded by the CMS detector in pp collisions at a 
centre-of-mass energy of 7 TeV. The sample corresponds to a total integrated luminosity of 
314 ± 34 nb _1 |T4| . During this data taking period, there were 1.6 pp collisions per bunch 
crossing, on average. ]/xp mesons are reconstructed in the decay channel. The event 

selection requires good quality data from the tracking, muon, and luminosity detectors, in 
addition to good trigger conditions. 

The analysis is based on events triggered by a double-muon trigger that requires the detection 
of two independent muon segments at LI, without any further processing at the HLT. All three 
muon systems, DT, CSC and RPC, take part in the trigger decision. The coincidence of two 
muon signals, without any cut on pj, is enough to keep the trigger rate reasonably low at the 
instantaneous luminosities of the LHC start-up. 

Events not coming from pp collisions, such as those from beam-gas interactions or beam- 
scraping in the transport system near the interaction point, which produce a large activity in 
the pixel detector, are removed by requiring a good primary vertex to be reconstructed |(T5| . 

3.2 Monte Carlo simulation 

Simulated events are used to tune the selection criteria, to check the agreement with data, 
to compute the acceptance, and to derive corrections to the efficiencies (Section [4b. Prompt 
J/tp mesons have been simulated using Pythia 6.421 [16], which generates events based on the 
leading-order color-singlet and color-octet mechanisms, with non-relativistic QCD (NRQCD) 
matrix elements tuned by comparing calculations with CDF data ||3] [17|. Color-octet states 
undergo a shower evolution. Simulated events with b-hadron pairs were also generated with 
Pythia and the b hadrons decayed inclusively into ]/ip using the EvtGen package [TBI . Final- 
state bremsstrahlung was implemented using PHOTOS [ 19112011 . 

The generated events were passed through the GEANT4-based [21] detector simulation and 
processed with the same reconstruction program as used for collision events. The detector 
simulation includes the trigger, as well as the effects of the finite precision of alignment and 
calibration, as determined using LHC collision data and cosmic-ray muon events l[22 | . 

3.3 Offline muon reconstruction 

In this analysis, muon candidates are defined as tracks reconstructed in the silicon tracker 
which are associated with a compatible signal in the muon chambers. 
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Two different muon reconstruction algorithms are considered [23 J. The first one provides high- 
quality and high-purity muon reconstruction for tracks with pj > 4 GeV/c in the central pseu- 
dorapidity region (\n \ < 1.3) and pj > 1 GeV/c in the forward region; these muons are referred 
to as Global Mnons. The second muon reconstruction algorithm achieves a better reconstruction 
efficiency at lower momenta; these muons are referred to as Tracker Muons. There is an overlap 
between these two reconstruction methods. If a muon is reconstructed by both algorithms, it 
is assigned to the Global Muon category alone, making the two categories exclusive. Global 
Muons have a higher reconstruction purity. In both cases, the track momentum is determined 
by the fit in the silicon tracker. 

To reduce muon backgrounds, mostly from decays in flight of kaons and pions, and to en- 
sure good quality reconstructed tracks, muon tracks are required to pass the following require- 
ments: they must have at least 12 hits in the tracker, at least two of which are required to be in 
the pixel layers, a track fit with a x 2 per degree of freedom smaller than four, and must pass 
within a cylinder of radius 3 cm and length 30 cm centered at the primary vertex and parallel to 
the beam line. If two (or more) tracks are close to each other, it is possible that the same muon 
segment or set of segments is associated with more than one track. In this case the best track is 
selected based on the matching between the extrapolated track and the segments in the muon 
detectors. 

The momentum measurement of charged tracks in the CMS detector has systematic uncer- 
tainties due to imperfect knowledge of the magnetic field, modeling of the detector material, 
sub-detector misalignment, and biases in the algorithms which fit the track trajectory; these 
effects can shift and /or broaden the reconstructed peaks of dimuon resonances. In addition 
to calibrations already applied to the data l|22l l24l |25| , residual effects can be determined by 
studying the dependence of the reconstructed dimuon peak shapes on the muon kinematics. 
The transverse momentum corrected for the residual scale distortion is parametrized as 

p™ = (l + a 1 +a 2 t ] 2 )pr as , (1) 

where p™ eas is the measured muon transverse momentum. A likelihood fit ||26| was performed 
to the invariant mass shapes to minimize the difference between the reconstructed ]/ip mass 
and the world-average value [27]. The resulting values of a\ and a-i are (3.8 ± 1.9) • 10~ 4 and 
(3.0 ± 0.7) • 10~ 4 , respectively. 

3.4 }lxp event selection 

To select the events with J/i/? decays, muons with opposite charge are paired and their invariant 
mass is computed. The invariant mass of the muon pair is required to be between 2.6 and 

3.5 GeV/ c 2 . The two muon trajectories are fitted with a common vertex constraint, and events 
are retained if the fit x 2 probability is larger than 0.1%. This analysis uses combinations of 
two Global Muons, two Tracker Muons, and one Global and one Tracker Muon. On average, 
1.07 J/tp combinations were found per selected dimuon event. In case of multiple combinations 
in the same event, the one with the purest muon content is chosen. If there are two or more 
dimuon candidates of the same type (Global-Global, Global-Tracker, or Tracker-Tracker) the 
one of highest pj is chosen. 

The opposite-sign dimuon mass spectrum is shown in Fig. [ljfor three different J/i/> rapidity 
ranges. About 27000 candidates have been reconstructed, of which about 19% are in the 
two-Global-Muon category, 54% in the Global-Tracker-Muon category, and the remaining in 
the two-Tracker-Muon category. 
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Figure 1: Opposite-sign dimuon invariant mass distributions in three ]/tp rapidity ranges, fitted 
with a Crystal Ball function plus an exponential (Section|5|. The poorer dimuon mass resolution 
at forward rapidity is caused by the smaller lever arm of the muon tracks. 

4 Acceptance and Efficiency 
4.1 Acceptance 

The acceptance reflects the finite geometrical coverage of the CMS detector and the limited 
kinematical reach of the muon trigger and reconstruction systems, constrained by the thickness 
of the material in front of the muon detectors and by the track curvature in the magnetic field. 



The ]/ip acceptance A is defined as the fraction of detectable ]/ip 
of the dimuon transverse momentum px arid rapidity y, 



}i + }i decays, as a function 



A(p T ,y; A B ) 



(2) 



N det(PT,y;A e ) 

N g en(pT,i/;A ) 

where N^ ei is the number of detectable ]/tp events in a given (pj, y) bin, expressed in terms 
of the dimuon variables after detector smearing, and Ng e n is the corresponding total number 
of generated ]/tp events in the Monte Carlo simulation. The parameter reflects the fact that 
the acceptance is computed for various polarization scenarios, as explained below. The large 
number of simulated events available allows the use of a much smaller bin size for determining 
A than what is used for the cross-section measurement. 

The criteria for detecting the muons coming from the J /xp decay is that both muons should be 
within the geometrical acceptance of the muon detectors and have enough momentum to reach 
the muon stations. The following kinematic cuts, defining the acceptance region, are chosen so 
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as to guarantee a single-muon detection probability exceeding about 10%: 

p£ > 3.3 GeV/c for |^| < 1.3 ; 

p ?< > 2.9 GeV/c for 1.3 < |/| < 2.2 ; 

p\ > 2.4 GeV/c for 2.2 < \rjV\ < 2.4 . 

To compute the acceptance, ]/tp events are generated with no cut on pj and within a rapidity 
region extending beyond the muon detector's coverage. 

The acceptance as a function of pj and \y\ is shown in the left plot of Fig. [2] for the combined 
prompt and non-prompt ]/xp mesons, with the prompt component decaying isotropically, cor- 
responding to unpolarized production. The right plot of Fig. [2] displays the pj and |y| distri- 
bution of muon pairs measured with an invariant mass within ± 100 MeV/ c 2 of the known 
J/ip mass 123. 




Figure 2: Left: Acceptance as a function of the J/tp pj and rapidity. Right: Number of muon 
pairs within ± 100 MeV/ c 2 of the nominal ]/ip mass, in bins of pi and \y\. 

Systematic uncertainties on the acceptance have been investigated, as described in the follow- 
ing paragraphs. 

• Final-state radiation. At the generator level, the dimuon momentum may differ 
from the J/xp momentum, due to final-state radiation (FSR). The difference between 
the acceptance computed using the dimuon system or the ]/tp variables in Eq. [2] is 
taken as a systematic uncertainty. 

• Kinematical distributions. Different spectra of the generated ]/xp might produce 
different acceptances. The difference between using the Pythia spectra and other 
theoretical calculations (mentioned in Section[7|) is taken as a systematic uncertainty. 

• b-hadron fraction and polarization. The J/t/> mesons produced in b-hadron decays 
can, in principle, have a different acceptance with respect to the prompt ones, due 
to their different momentum spectra, leading to an uncertainty coming from the un- 
known proportion of b hadrons in the inclusive sample. This fraction has been var- 
ied in the Monte Carlo simulation by 20%, the average accuracy of the measurement 
performed here (presented in Section|6]l; the difference between the two acceptances 
is taken as an estimate of this uncertainty. For non-prompt ]/ip mesons the b-hadron 
events are generated with the ]/ip polarization as measured by the BaBar experi- 
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4 Acceptance and Efficiency 



ment |28| , and the corresponding systematic uncertainty is evaluated by taking the 
difference with respect to the one predicted by EvtGen. 

• px calibration and resolution. A difference between the muon momentum scale in 
data and simulated events would lead to a different acceptance. The muon trans- 
verse momenta have been calibrated as described in Section 13.31 The maximum 
residual bias remaining after the calibration is estimated to be 0.05%. As a con- 
servative estimate, a bias equivalent to this residual uncertainty is applied to the 
simulated muon momenta. The change in the recomputed acceptance is taken as a 
systematic uncertainty. Similarly a difference in the momentum resolution between 
data and simulated events would also give a different acceptance. The acceptance 
has been computed with simulated muon momenta smeared according to the reso- 
lution measured in data [26 J and the difference is taken as a systematic uncertainty. 

Finally the distribution of the z position of the pp interaction point could in principle influence 
the acceptance. Several Monte Carlo samples of }/ip mesons have been generated, each coming 
from different positions along the beam line (between —10 and +10 cm with respect to the 
centre of the collision region) and a negligible variation of the acceptance has been found. 

4.2 Efficiency 

The single-muon efficiency is computed using the Tag-and-Probe method |23ll29| . The combined 
trigger and offline-reconstruction efficiency for a single muon is measured with a data sample 
collected with looser trigger requirements and is defined as 

£ W = e track ' £ id | track ' e trig | track+id ' ( 3 ) 

where e^rack * s * ne tracking efficiency, | trac ] c is the muon identification efficiency in the 
muon systems for a tracker-reconstructed muon, and finally e tri g | track+id * s * ne P roDa b m ty 
for an offline reconstructed muon to have also fired the trigger. 

The tracking efficiency is constant in the momentum range defined by the acceptance cuts, and 
it varies only slightly in the <p — t] plane [29]. The muon identification and trigger efficiencies 
have a stronger pj and \r]P\ dependence, which is mapped with a finer granularity (nine to 
twelve pj and five \rj^\ bins). 

The efficiency to detect a given J/xp event is thus dependent on the value of the muon-pair 
kinematic variables, and is given by 

e(J/V) = ■ •(! + *>)■ ^vertex ■ (4) 



The factor p represents a correction to the factorization hypothesis and is evaluated from the 
Monte Carlo simulation. The non-vanishing values of p, varying between —0.19 and 0.30, are 
mainly due to the relatively large bin sizes used to determine the muon efficiencies. 

The efficiency for the two muon tracks to be consistent with coming from a common vertex 



(Section 3.4 1, % e rtex/ is measured to be (98.35 ± 0.16)%, by comparing the number of two- 
Global-Muon combinations within ±100 MeV/c 2 of the nominal ]/xp mass with and without 
the common vertex requirement. Given the precision of this estimate, the corresponding sys- 
tematic uncertainty can be neglected. The following systematic uncertainties on the J/tp effi- 
ciency are considered: 

• p factor. Any variation of the muon spectrum within each large bin may lead to a 
different value of p. By reweighting the Pythia Monte Carlo simulation, we vary the 
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]/tp pj spectrum to reproduce different theoretical predictions (Section [7), and take 
the largest variation as the systematic uncertainty on p. 

• Muon efficiency. The statistical uncertainty on each muon efficiency is propagated 
using toy Monte Carlo experiments, and the r.m.s. of the newly computed ]/tp effi- 
ciencies are assigned as systematic uncertainties. The largest systematic errors are in 
the bins with less events or in those where the background is largest. When selecting 
the tag muon, the Tag-and-Probe method produces a slight bias on the kinematics 
of the probe muon, hence a small difference arises between the measured single- 
muon efficiencies and those of an unbiased sample. This small effect is studied in 
the Monte Carlo simulation and corrected for. The whole correction is conservatively 
taken as a systematic uncertainty on the efficiencies and summed in quadrature with 
the statistical uncertainty. 

5 Inclusive ]lxp cross section 

The measurement of the inclusive px differential cross section is based on the equation 



dp T dy r r r r j J Ldt ■ Ap T ■ Ay 

where NcorrQ/VO i s the ]/xp yield, corrected for the ]/ip acceptance and selection efficiency, in 
a given transverse momentum-rapidity bin, J Ldt is the integrated luminosity, Apx and Ay 
are the sizes of the px and rapidity bins, and BRQ/xp — > ]d + fi~) is the branching ratio of the 
]/\p decay into two muons. 

5.1 }lxp yields 

The corrected yield, NcorrQ/^), is determined in two steps. First, in each rapidity and px bin 
an unbinned maximum likelihood fit to the invariant mass distribution is performed. 

The resulting yield is then corrected by a factor that takes into account the average acceptance 
(A) and detection efficiency (e) in the bin under consideration. 

In the mass fits, the shape assumed for the signal is a Crystal Ball function [30J, which takes into 
account the detector resolution as well as the radiative tail from bremsstrahlung. The shape of 
the underlying continuum is described by an exponential. Table [l] lists the J/i/? uncorrected 
signal yields and the corresponding statistical uncertainties from the fit, for the chosen bins. 

Different functions were used to assess systematic effects coming from the fit function chosen to 
model the signal and the continuum shapes. For the signal, the Crystal Ball function was varied 
to a sum of a Crystal Ball and a Gaussian, while for the background a second-order polynomial 
was used. The maximum difference in the result was taken as a systematic uncertainty. The 
uncertainty is particularly large for the low-px bins, where the signal purity is the smallest. 

Additionally, a bias on the muon momentum scale can shift the events from one ]/ip px bin 
to the adjacent ones. To estimate this systematic effect, a bias has been applied to the muon 
momenta equal to the residual uncertainty on the scale after the calibration, as explained in 
Section [3^4} and a negligible variation was found. 

5.2 Inclusive }lxp cross section results 

The previously discussed systematic uncertainties affecting the inclusive ]/ip cross section are 
listed in Table [2] In addition, the relative error on the luminosity determination is 11%, and is 
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5 Inclusive Jfy cross section 



Table 1: Uncorrected event yield (with its statistical error from the fit) in each pj bin, together with the 
average acceptance times efficiency (computed in the unpolarized production scenario). 
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\y\ < i-2 


0.50- 


-0.75 


829.3 ±44.7 


0.079 ± 0.010 


6.5- 


8.0 


726.5 ± 28.3 


0.084 ± 0.005 


0.75- 


-1.00 


1006.0 ± 48.8 


0.078 ± 0.010 


8.0- 


10.0 


868.1 ±30.7 


0.178 ± 0.005 


1.00- 


- 1.25 


1216.8 ±52.8 
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0.075 ± 0.008 
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-2.00 


1132.7 ±57.5 
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-2.25 


1122.7 ± 55.0 


0.071 ± 0.006 


1.2 < \y\ < 1.6 
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903.3 ± 72.4 
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401.7 ±23.2 


0.035 ± 0.004 


2.75- 


-3.00 


757.6 ± 36.2 


0.077 ±0.005 


4.5- 


5.5 


618.9 ± 28.9 


0.086 ± 0.004 
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-3.25 


756.1 ± 35.7 


0.082 ± 0.005 


5.5- 


6.5 


690.9 ± 34.0 


0.167 ± 0.005 


3.25- 


-3.50 


703.6 ± 33.6 


0.084 ± 0.004 


6.5- 


8.0 


712.0 ± 28.0 


0.247 ± 0.006 


3.50- 


-4.00 


1150.2 ±40.0 


0.092 ± 0.005 


8.0- 


10.0 


463.7 ±23.3 


0.334 ± 0.009 


4.00- 


-4.50 


991.8 ± 35.8 


0.100 ± 0.004 


10.0- 
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406.2 ± 22.4 


0.445 ± 0.010 
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1441.4 ±42.6 


0.117 ±0.005 
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Table 2: Relative systematic uncertainties on the corrected yield for different J/ip rapidity bins. The 
variation range over the different pj bins is given. In general, uncertainties depend only weakly on the 
pi values, except for the fit function systematic uncertainty, which decreases with increasing pj due to 
the better purity of the signal. The large excursion of the muon efficiency systematic uncertainty reflects 
changes in the event yield and in the signal purity among the pj bins. 
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Relative error (%) 










\y\ 


< 1.2 


1.2 < \y\ < 1.6 1.6 < 


y\ < 2.4 


Acceptance 


FSR 


0.8 


-2.5 


0.3 - 1.6 


0.0 


-0.9 




pj calibration and resolution 


1.0 


-2.5 


0.8 - 1.2 


0.1 


-1.0 




Kinematical distributions 


0.3 


-0.8 


0.6 - 2.6 


0.9 


-3.1 




b-hadron fraction and polarization 


1.9 


-3.1 


0.5 - 1.2 


0.2 


-3.0 


Efficiency 


Muon efficiency 


1.9 


-5.1 


2.3 - 12.2 


2.7 


-9.2 




p factor 


0.5 


-0.9 


0.6 - 8.1 


0.2 


-7.1 


Yields 


Fit function 


0.6 


- 1.1 


0.4 - 5.3 


0.3 


-8.8 



common to all bins. Table [3] reports the values of the resulting ]/\p differential cross section, for 
different polarization scenarios: unpolarized, full longitudinal polarization and full transverse 
polarization in the Collins-Soper or the helicity frames f7\ . 

Figure |ij shows the inclusive differential cross section d d ^ dy ■ BR(]/xp —> in the three 

rapidity ranges, showing statistical and systematic uncertainties, except the luminosity uncer- 
tainty, added in quadrature. It should be noted that the first bin in the forward rapidity region 
extends down to zero ]/ip pj. 

The total cross section for inclusive ]/ip production, obtained by integrating over pj between 
6.5 and 30GeV/c and over rapidity between —2.4 and 2.4, in the unpolarized production hy- 



5.2 Inclusive J/xp cross section results 
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6 Fraction of J/i/> from b-hadron decays 
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Figure 3: Differential inclusive J/i/? cross section as a function of px for the three different ra- 
pidity intervals and in the unpolarized production scenario. The errors on the ordinate values 
are the statistical and systematic errors added in quadrature. The 11% uncertainty due to the 
luminosity determination is not shown and is common to all bins. 

pothesis, gives 

cr(pp -> ]/lp + X) • BR(J/t/> -> y+jT ) = 97.5 ± 1.5(stat) ± 3.4(syst) ± 10.7(luminosity) nb. (6) 



6 Fraction of J/tp from b-hadron decays 

The measurement of the fraction of ]/xp yield coming from b-hadron decays relies on the dis- 
crimination of the J/tp mesons produced away from the pp collision vertex, determined by the 
distance between the dimuon vertex and the primary vertex in the plane orthogonal to the 
beam line. 

The primary vertices in the event are found by performing a common fit to tracks for which the 
points of closest approach to the beam axis are clustered in z, excluding the two muons forming 
the ]/tp candidate and using adaptive weights to avoid biases from displaced secondary ver- 
tices. Given the presence of pile-up, the primary vertex in the event is not unique. According 
to Monte Carlo simulation studies, the best assignment of the primary vertex is achieved by 
selecting the one closest in the z coordinate to the dimuon vertex. 



6.1 Separating prompt and non-prompt J/ip 
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6.1 Separating prompt and non-prompt J/tp 

As an estimate of the b-hadron proper decay length, the quantity £yy = L xy ■ myy/ pj is com- 
puted for each ]/tp candidate, where my^ is the ]/tp mass f27\ and L xy is the most probable 
transverse decay length in the laboratory frame [|3Tll32| . L xy is defined as 

_ u T a he 

where x is the vector joining the vertex of the two muons and the primary vertex of the event, 
in the transverse plane, u is the unit vector of the J/t/> pj, and a is the sum of the primary and 
secondary vertex covariance matrices. 

To determine the fraction /g of J /xp mesons from b-hadron decays in the data, we perform an 
unbinned maximum-likelihood fit in each pj and rapidity bin. The dimuon mass spectrum and 
the £}/y distribution are simultaneously fit by a log-likelihood function, 



N 

lnL = X>f(%,m w ) , (8) 

i=i 



where N is the total number of events and m«„ is the invariant mass of the muon pair. The 
expression for F(£yy>, m y}l ) is 

F(%, m w ) = fsi g ■ F S i g {e ]/t p) ■ M Sis (m m ) + (1 - f Sig ) • F B fc g (%) ■ M Bkg (m m ) , (9) 

where: 

• fsig is the fraction of events attributed to ]/xp sources coming from both prompt and 
non-prompt components; 

• Msi g (m}i}i) and Mb^wj^) are functional forms describing the invariant dimuon 
mass distributions for the signal and background, respectively, as detailed in Sec- 
tion [51} 

• Fsi g (£]/y) and FBicg(^j/i/j) are functional forms describing the £yy distribution for the 
signal and background, respectively. 

The signal part is given by a sum of prompt and non-prompt components, 

Fsi g (£ J / f )=fB-F B (£y jp ) + (l-f B )-F p (£y ip ) , (10) 

where fs is the fraction of // ip from b-hadron decays, and Fp(£yy) and Fb(%i/0 are 
the £yy distributions for prompt and non-prompt // ip, respectively. 

As £yy should be zero in an ideal detector for prompt events, F p (£yy) is described 
simply by a resolution function. The core of the resolution function is taken to be a 
double-Gaussian and its parameters are allowed to float in the nominal fit. Since £yy 
depends on the position of the primary vertex, an additional Gaussian component 
is added, to take into account possible wrong assignments of the primary vertex; its 
parameters are fixed from the Monte Carlo simulation. 



The £yy shape of the non-prompt component in Eq. 10 is given by convolving the 
same resolution function with the true £yy distribution of the ]/xp from long-lived b 
hadrons, as given by the Monte Carlo simulation. 
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6 Fraction of J/i/> from b-hadron decays 



For the background tyy distribution FskgiZyq), the functional form employed by 
CDF is used: 

F Bkg (x) = (l-/ + -/_-/ sym )R(*) + 

A+ A_ 



/sym 
2A S y m 



e As y m ] <S>R(x' - x) , (11) 



where R(x) is the resolution model mentioned above, (z = {+, — ,sym}) are the 
fractions of the three long-lived components with mean decay lengths A,-, and 6{x) 
is the step function. The effective parameters A, are previously determined with a fit 
to the £-j/ip distribution in the sidebands of the dimuon invariant mass distribution, 
defined as the regions 2.6-2.9 and 3.3-3.5 GeV/c 2 . 

The parameter /g (b fraction) is determined in the same rapidity regions as used to present the 
inclusive production cross section but some pj bins are grouped, since more events per bin are 
needed to determine all fit parameters. Figure [4] shows the projection of the likelihood fits in 
two sample bins. The full results are reported in Table |3J where fs has been corrected by the 
prompt/ non-prompt acceptances, as discussed in Section |4j The fitting procedure has been 
tested in five sample bins using toy experiments, which establish reasonable goodness-of-fit 
and exclude the possibility of biases in the /# determination. 




Figure 4: Projection in the £jm dimension of the two-dimensional likelihood fit (in mass and lyy) 
in the bins 2 < p T < 4.5GeV/c, 1.2 < \y\ < 1.6 (left) and 6.5 < p T < lOGeV/c, 1.6 < \y\ < 2.4 
(right), with their pull distributions (bottom). 

Figure [5] shows the measured b fraction. It increases strongly with pj. At low pj, essentially all 
]/tp mesons are promptly produced, whereas at pj ~ 12GeV/c around one third come from 



6.1 Separating prompt and non-prompt J/ip 
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Table 4: Fit results for the determination of the fraction of J/ip mesons from b hadrons in pj and 
|y| bins, corrected by the prompt and non-prompt acceptances. The average pj and r.m.s. per 
bin are also quoted. The two uncertainties in the b-fraction values are statistical and systematic, 
respectively. 
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beauty decays. This pattern does not show a significant change with rapidity (within the 
current uncertainties) over the window covered by the CMS detector. The CMS results are 
compared to the higher-precision data of CDF [5], obtained in proton-antiproton collisions at 
^/s = 1.96 TeV. It is interesting to note that the increase with pj of the b fraction is very similar 
between the two experiments, the CMS points being only slightly higher, despite the different 
collision energies. 
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Figure 5: Fraction of the ]/\p production cross section originating from b-hadron decays, as 
a function of the }/\p pj, as measured by CMS in three rapidity bins and by CDF, at a lower 
collision energy. 
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6 Fraction of J/i/> from b-hadron decays 



6.1 .1 Systematic uncertainties affecting the b-f raction result 

Several sources of systematic uncertainty have been addressed and are described in the follow- 
ing lines. 

• Residual misalignments in the tracker. The effect of uncertainties in the measured 
misalignment of the tracker modules is estimated by reconstructing the data several 
times using different sets of alignment constants. These sets reflect the uncertainty in 
the constants and, in particular, explore possible deformations of the tracker which 
are poorly constrained by the data Il22| . The largest difference between the results 
with the nominal set of constants and with these sets is taken as a systematic uncer- 
tainty. 

• b-hadron lifetime model. In an alternative approach, £yy is described by a convolu- 
tion of an exponential decay with a Gaussian function, which describes the smearing 
due to the relative motion of the J/t/> with respect to the parent b hadron. The differ- 
ence between the nominal Monte Carlo template model and this alternative is taken 
as a systematic uncertainty. 

• Primary vertex estimation. In an alternative approach, the beam spot as calculated 
on a run-by-run basis is chosen as the primary vertex in calculating £j/p and the fit 
is repeated. The difference is taken as a systematic uncertainty. 

• Background. The background is fitted using only the sidebands and the result is 
used as input to the fit in the signal region. The effect of a ± 100 MeV/ c 2 variation 
in the sideband boundaries is taken as a systematic uncertainty. 

• lyy resolution model. The nominal (triple-Gaussian) fit model for the decay length 
resolution is compared to a model using two Gaussians only, fixing the "additional" 
Gaussian to be zero. The difference is taken as a systematic uncertainty. 

• Different prompt and non-prompt efficiencies. The Monte Carlo simulation pre- 
dicts small differences between the prompt and non-prompt ]/xp efficiencies. These 
are taken into account and the relative difference assumed as a systematic uncer- 
tainty. 

A summary of all systematic effects and their importance is given in Table [5] 

Table 5: Summary of relative systematic uncertainties in the b-fraction yield (in %). The vari- 
ation range over the different pj bins is given in the three rapidity regions. In general, uncer- 
tainties are px~dependent and decrease with increasing pj. 



\y\<12 1.2<|y|<1.6 1.6 < |y| < 2.4 



Tracker misalignment 


0.5 


-0.7 


0.9- 


-4.6 


0.7- 


-9.1 


b-lifetime model 


0.0 


-0.1 


0.5- 


-4.8 


0.5- 


- 11.2 


Vertex estimation 


0.3 


1.0- 


- 12.3 


0.9- 


-65.8 


Background fit 


0.1 


-4.7 


0.5- 


-9.5 


0.2- 


- 14.8 


Resolution model 


0.8 


-2.8 


1.3- 


- 13.0 


0.4- 


-30.2 


Efficiency 


0.1 


- 1.1 


0.3- 


-1.3 


0.2- 


-2.4 



6.1 .2 Prompt and non-prompt J/tp production cross sections 

The prompt ]/ip cross section and the cross section from b-hadron decays, together with their 
statistical and systematic uncertainties, are given in Tables |6]and|7j respectively, for the different 
polarization scenarios considered in Section |5j 

The total cross section for prompt ]/ip production times BR(]/ip — > for the unpolarized 
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Table 6: Differential prompt J/tp cross sections for each polarization scenario considered: un- 
polarized (Ag = 0), full longitudinal polarization (Ag = —1) and full transverse polarization 
(Ag = +1) in the Collins-Soper (CS) or the Helicity (HX) frames [7\. For the unpolarized case, 
the first error is statistical and the second is systematic; for the others the total error is given. 



Pt 
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0.134 ±0.033 ±0.016 
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149.1 ± 22.0 


86.7 ±12.8 


140.7 ± 20.8 


2.00- 


-2.75 


93.7 ±3.4 ±12.4 


65.8 ±9.1 


118.8 ± 16.3 


72.7 ±10.0 


110.3 ±15.2 


2.75- 


-3.50 


62.6 ±2.0 ±7.9 


44.5 ± 5.7 


78.8 ± 10.2 


49.1 ±6.4 


72.7 ±9.5 


3.50- 


-4.50 


37.4 ±1.1 ±4.9 


27.4 ± 3.7 


45.7 ± 6.2 


29.9 ±4.1 


42.8 ± 5.8 


4.50- 


-6.50 


15.2 ±0.4 ±2.0 


11.9 ± 1.6 


18.0 ± 2.4 


12.6 ±1.7 


17.1 ±2.3 


6.50- 


- 10.00 


3.08 ±0.11 ±0.37 


2.79 ± 0.35 


3.36 ± 0.42 


2.64 ± 0.33 


3.37 ± 0.42 


10.00 


- 30.00 


0.093 ±0.007 ±0.012 


0.092 ±0.014 


0.096 ± 0.014 


0.082 ±0.012 


0.100 ±0.015 



production scenario, has been obtained by integrating the differential cross section over pj 
between 6.5 and 30 GeV/c and over rapidity between —2.4 and 2.4, 

BR(]/ip -> p+ft - ) '^(PP -> prompt J/ip) = 70.9 ±2.1 ±3.0 ±7.8 nb , (12) 

where the three uncertainties are statistical, systematic and due to the measurement of the 
integrated luminosity, respectively. Similarly, the cross section of non-prompt ]/xp mesons from 
b-hadron decays, times BR(]/xp — > is 

BR(]/xp^p + ji-)-a(pp^bX^]/xpX) =26.0 ±1.4 ±1.6 ±2.9 nb . (13) 

The sum of these two cross sections differs slightly from the inclusive value, which was deter- 
mined assuming a b fraction taken from Monte Carlo expectations. 

7 Comparison with theoretical calculations 

The prompt ]/xp differential production cross sections, in the rapidity ranges considered in the 
analysis, as summarized in Table|6j were compared with calculations made with the Pythia |T6| 
and CASCADE J33J|34| event generators, as well as with the Color Evaporation Model (CEM) [35- 
39 1 . These calculations include the contributions to the prompt ]/tp yield due to feed-down 
decays from heavier charmonium states (x c and xp(2S)) and can, therefore, be directly com- 
pared to the measured data points, as shown in Fig. |6j In contrast, it is not possible to compare 
our measurement with the predictions of models such as the Color-Singlet Model (including 
higher-order corrections) I40j-l43| or the LO NRQCD model (which includes singlet and octet 
components) [44, 45 j, because they are only available for the direct }/ip production component, 
while the measurements include a significant contribution from feed-down decays, of the order 
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Table 7: Differential non-prompt J/t/> cross section times the ]/ip branching ratio to dimuons, 
assuming the polarization measured by the BaBar experiment [28J at the Y(4S). The first uncer- 
tainty is statistical and the second is systematic. 



p^ BKCJ/V-^-)-- 
(GeV/c) (nb/GeV/c) 



\v\ < i-2 



6.5- 


10.0 


1.30 ±0.08 ±0.19 


10.0- 


-30.0 


0.087 ±0.024 ±0.010 


1.2 < \y\ < 1.6 


2.0- 


4.5 


8.67 ±1.36 ±2.71 


4.5- 


6.5 


4.04 ± 0.41 ± 0.79 


6.5- 


10.0 


0.98 ±0.09 ±0.11 


10.0- 


-30.0 


0.065 ± 0.007 ± 0.008 


1.6 < \y\ < 2.4 


0.00- 


- 1.25 


4.31 ± 1.59 ± 3.54 


1.25- 


-2.00 


11.0 ±1.8 ±4.2 


2.00- 


-2.75 


11.9 ±1.4 ±3.4 


2.75- 


-3.50 


10.1 ±1.1 ±1.6 


3.50- 


-4.50 


7.19 ± 0.65 ± 1.25 


4.50- 


-6.50 


3.28 ± 0.24 ± 0.53 


6.50- 


- 10.00 


0.95 ±0.07 ±0.13 


10.00 


- 30.00 


0.055 ± 0.005 ± 0.007 



of 30% Il46ll47) . At forward rapidity and low pj the calculations underestimate the measured 
yield. 

The non-prompt ]/xp differential production cross sections, as summarized in Table[7j have been 
compared with calculations made with the Pythia and CASCADE Monte Carlo generators, and 
in the FONLL framework flQ|. The measured results are presented in Fig. [7] and show a good 
agreement with the calculations. 

8 Conclusions 

We have presented the first measurement of the J/tp production cross section in pp collisions 
at = 7 TeV, based on 314 nb _1 of integrated luminosity collected by the CMS experiment 
during the first months of LHC operation. 

The px differential J/t/> production cross section, in the dimuon decay channel, has been mea- 
sured in three rapidity ranges, starting at zero px for 1.6 < |y| < 2.4, at 2 GeV/c for 1.2 < |y| < 
1.6, and at 6.5 GeV/c for |y| < 1.2. The measured total cross section for prompt ]/ip production 
in the unpolarized scenario, in the dimuon decay channel, is 

a{pp -> ]/xp + X) ■ BR(J/i/> -> p + p~ ) = 70.9 ± 2.1(stat) ± 3.0(syst) ± 7.8(luminosity) nb , 

for transverse momenta between 6.5 and 30GeV/c and in the rapidity range \y\ < 2.4. Aside 
from the luminosity contribution, the systematic uncertainty is dominated by the statistical 
precision of the muon efficiency determination from data. 

The measured total cross section times BR(]/tp — > for ]/tp production due to b-hadron 



17 



CD 

CD 

JO 



10 



% 1 
O 

C\J 

"O 



10' 



I I 


— CMS data ] 
PYTHIA 




— CASCADE j 
CEM 




|y |<1.2 1 


r CMS,\/s = 7TeV 




: L = 314nb" 1 

i . . i , , , , i 





— CMS data 
PYTHIA 

— CASCADE 
CEM 



10 



15 20 25 30 

n J/>t 




pf (GeV/c) 



15 20 25 30 

J/V / 



p; ¥ (GeV/c) 




15 20 25 30 
p^" (GeV/c) 

Figure 6: Differential prompt J/i/> production cross section, as a function of px f° r the three different 
rapidity intervals. The data points are compared with three different models, using the PYTHIA curve 
to calculate the abscissa where they are plotted 1481 . 

decays, for 6.5 < pj < 30 GeV/c and \y\ < 2.4, is 

cr(pp ->bX-> ]/xpX) ■ BRQ/xp -> = 26.0 ± 1.4 (stat) ± 1.6 (syst) ± 2.9 (luminosity) nb . 



The differential prompt and non-prompt measurements have been compared with theoretical 
calculations. A reasonable agreement is found between data and theory for the non-prompt 
case while the measured prompt ]/ip cross section exceeds the expectations at forward rapidity 
and low pj. 
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